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Abstract.  Subsurface  optical  layers  distributed  at  two  different  depths  were  investigated  in 
Monterrey  Bay,  East  Sound,  and  the  Black  Sea  based  on  spatial  statistics  of  remote  sensing 
reflectance  (/?,*).  The  main  objective  of  this  study  was  to  evaluate  the  use  of  /?rs(443)//?rs(490) 
(hereafter  R 1)  skewness  (VO  as  an  indicator  of  vertical  optical  structure  in  different  marine 
regions.  Measurements  of  inherent  optical  properties  were  obtained  using  a  remotely  op¬ 
erated  towed  vehicle  and  R\  was  theoretically  derived  from  optical  profiles.  Although  the 
broad  range  of  trophic  status  and  water  stratification,  a  common  statistical  pattern  consist¬ 
ing  of  lower  xf/R  1 — a  deeper  optical  layer  was  found  in  all  study  cases.  This  variation  was 
attributed  to  optical  changes  above  the  opticline  and  related  to  horizontal  variability  of  par¬ 
ticulates  and  spectral  variations  with  depth.  We  recommend  more  comparisons  in  stratified 
coastal  waters  with  different  phytoplankton  communities  before  the  use  of  \//R  1  can  be  gener¬ 
alized  as  a  noninvasive  optical  proxy  for  screening  depth  changes  on  subsurface  optical  layers. 
© 2011  Society  of  Photo-Optical  Instrumentation  Engineers  (SPIE).  [DOI:  10. 1 1 17/1.3634055] 
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1  Introduction 

A  common  assumption  of  remote  sensing  algorithms  based  on  ocean  color  sensors  is  the  vertical 
homogeneity  of  the  water  column  in  terms  of  optical  properties.  This  approximation  is  very  often 
not  met  in  coastal  and  oceanic  stratified  waters  due  to  the  presence  of  laminar  features  altering 
the  underwater  light  field.  These  submerged  layers  commonly  correspond  with  thin  layers  (i.e., 
<3  m  thick),  have  typically  high  concentration  of  dissolved  and  particulate  material  with  respect 
to  the  surrounding  medium,  and  are  preferentially  developed  along  the  horizontal  component.1 
Technologies  used  to  detect  vertical  location  of  subsurface  optical  layers  are  generally  based  on 
lidar  (light  and  range  detection)  systems.1-'  Unlike  these  investigations,  in  a  recent  contribution 
we  showed  a  new  approach  to  discriminate  waters  with  shallow  versus  deep  optical  layers  based 
on  passive  optical  measurements.4  Briefly,  the  relative  distance  of  the  optical  layer  to  the  sea 
surface  is  estimated  by  calculating  the  third  moment  around  the  mean  (i.e.,  skewness  or  VO 
of  a  specific  remote  sensing  reflectance  ratio  [R\  =  /?rs(443)//?rs(490)].  As  x//R\  decreases  the 
subsurface  optical  layer  becomes  deeper. 
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The  present  study  has  four  objectives.  First,  it  will  test  optical  relationships  previously 
found  by  Ref.  4  by  using  an  alternative  approach  consisting  of  independent  calculations  of 
x//R\  based  on  in  situ  measurements  of  inherent  optical  properties  (lOPs)  obtained  with  an 
undulating  remotely  operated  towed  vehicle  (ROTV).  Second,  it  will  examine  whether  the 
aforementioned  approach  can  be  generalized  across  different  marine  coastal  domains  located 
at  different  latitudes  and  characterized  by  different  trophic  status.  Third,  it  will  investigate  the 
mechanisms  explaining  V^l  changes  in  terms  of  IOP  modifications.  Finally,  it  will  quantify 
the  influence  of  water  stratification  on  biological  “layering”  and  its  impact  on  \//R\  patterns. 

2  Methods 

2.1  ROTV  Surveys 

High  vertical  resolution  profiles  (0.02  to  0.56  m)  of  in  situ  optical  (i.e.,  total  absorp¬ 
tion  coefficient  of  dissolved  +  particulate  matter,  a,  and  beam  attenuation  coefficient, 
c)  and  conductivity-temperature-depth  (CTD)  measurements  were  obtained  using  a  ROTV 
(Scanfish  MK  II,  intelligent  undulating  vehicle,  EIVA,  Denmark)  in  coastal  oligotrophic 
[Black  Sea  (BS),  41.22°N  to  43.70°N,  28.90°E  to  31.26°E],  mesotrophic  [Monterrey  Bay 
(MB),  36.33°N  to  36.83°N,  121.05°W  to  1 22.85° W],  and  eutrophic  [East  Sound  (ES), 
48.62°N  to  48.67°N,  1 22.86° W  to  1 22.89° W)  waters  during  August  to  September  2008, 
June  2008,  and  May  2009,  respectively.  The  ROTV  can  provide  higher  temporal  (up  to 
1 0-fold)  and  spatial  (up  to  25-fold  vertical,  up  to  1 0-fold  horizontal)  resolved  cross  sections  than 
oceanographic  gliders  (e.g.,  electric  Slocum).  Thus,  and  given  our  future  interest  in  implement¬ 
ing  a  satellite-derived  approach  for  detecting  depth  of  optical  marine  layers,  we  chose  Scanfish 
as  the  sampling  tool  in  this  investigation. 

The  ROTV  (0.9  x  0.3  x  1 .8  m,  length,  height,  width)  was  operated  using  an  undulating 
mode  complemented  with  a  winch  to  expand  the  vertical  sampling  range  up  to  400  m.  The 
ROTV  has  a  weight  of  75  kg  and  a  payload  of  approximately  50  Kg.  ROTV  data  were  collected 
between  14  and  1 6  p.m.  (local  time)  and  two  case  studies  were  selected  per  study  site:  “shallow” 
(NS)  and  “deep”  (FS)  subsurface  optical  layer. 


2.2  Analysis  of  Optical  and  CTD  Profiles 

Raw  a  and  c  coefficients  were  derived  from  an  ac-9  (WetLabs)  at  three  wavelengths  (440,  488, 
and  675  nm)  and  corrected  by  salinity  and  temperature  effects5,6  using  CTD  (SeaBird  911) 
measurements.  Scattering  residuals  were  removed.7  Only  descending  dives  including  ROTV- 
derived  IOPs,  temperature,  and  conductivity  were  processed  and  smoothed  every  1  m  along 
the  vertical.  Salinity  and  seawater  density  at  each  depth  was  derived  from  CTD  variables  and 
using  the  standard  UNESCO  polynomial  equation  of  state. 8,g  For  each  variable,  data  gaps  were 
removed  by  linear  interpolation.  Missing  determinations  were  more  common  at  finer  spatial 
sampling  frequencies  (e.g.,  in  ES).  The  pycnocline  depth  was  defined  as  the  depth  at  which 
vertical  changes  of  water  density  were  above  0.01 5  Kg  m~4. 


2.3  Modeling  of  Ocean  Color 

Spectral  remote  sensing  reflectance  [/?rs(^)]  was  estimated  for  each  descending  profile  based 
on  vertical  distribution  of  a  and  c  values  at  specific  wavelengths.  A  light  propagation  model 
(Hydrolight,  Sequoia  Inc.)  and  ancillary  data  provided  from  airport-based  meteorological 
stations  (Black  Sea,  www.infospace.ru,  Monterrey  and  Eastsound,  ww  w.wunderground.com) 
were  used  to  simulate  /?rs(440)  and  RTS( 488)  above  each  Scanfish  profile.  For  each  case  study, 
we  computed  the  skewness  of  /?rs(440)/J?rs(488)  ratios4  and  its  uncertainty  by  creating  synthetic 
Scanfish  transects  with  the  same  horizontal  length  (i.e.,  n  =  10),  but  using  different  combinations 
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Fig.  1  An  example  of  Scanfish  MK  II  dives  (ups  and  downs)  in  the  BS,  MB,  and  ES.  BS,  MB  and 
ES  sampling  locations  had  an  averaged  bottom  depth  of  125,  120,  and  27  m,  respectively. 


of  dives  (ups  and  downs)  obtained  randomly  from  the  original  datasets.  For  each  experiment, 
we  examined  how  the  skewness  of  /?rs(440)//?rs(488)  ratios  and  optical  properties  measured 
above  the  subsurface  layer  [e.g.,  a(675)]  were  interrelated.  Also,  we  explored  mechanisms 
behind  these  statistical  relationships  by  analyzing  the  response  of  \//R\  to  magnitude  changes 
on  a  and  total  scattering  (i.e.,  b  =  c  —  a)  coefficients  at  different  wavelengths. 

Finally,  we  investigated  the  depth  of  the  signal  source  based  on  first  optical  depth  (FOD) 
estimates.  FOD  was  calculated  as  the  inverse  of  488)  or  the  diffuse  attenuation  coefficient 
of  downwelling  irradiance,  with  A^/(488)  ~  a(488)//7  (Ref.  10),  where  jl  is  the  average  cosine 
of  the  zenith  angle  of  refracted  solar  photons  (direct  beam)  just  beneath  the  sea  surface. 


3  Results  and  Discussion 

The  ROTV  diving  “saw  tooth”  pattern  varied  between  studied  areas  (Fig.  1).  Depth  range, 
averaged  vertical  resolution,  and  horizontal  spacing  between  profiles  were  2.8  to  65,  0.41,  and 
370.9  m  for  BS,  1.6  to  48,  0.36,  and  268.5  m  for  MB,  and  0.7  to  22,  0.02,  and  87.4  m  for  ES, 
respectively.  Cross-sections  of  a(615)  suggested  that  phytoplankton  was  an  important  optical 
component  in  all  subsurface  layers  under  investigation  (Fig.  2).  Consistent  with  an  increase 
of  solar  radiation  and  freshwater  river  discharge  as  the  spring-summer  season  progresses,  the 
pycnocline  and  the  optical  subsurface  layer  characterized  by  relatively  high  a{615)  with  respect 
to  background  (up  to  3-fold)  were  always  deeper  during  the  second  experiment  in  chronological 
order.  This  variation  was  remarkable  in  BS,  where  the  vertical  position  of  maximum  a(675)  and 
pycnocline  was  14  and  8.8  m  in  NS,  and  28  and  23.2  m  in  FS,  respectively.  In  general,  drastic 
vertical  changes  of  a(675)  were  observed  in  the  vicinity  of  the  pycnocline  but  during  the  FS 
experiment  in  ES  [Fig.  2(f)].  In  this  case,  it  is  suggested  that  phytoplankton  communities  were 
not  actively  growing  and  were  probably  sinking  as  part  of  a  post-bloom  stage  mainly  composed 
of  senescent  cells  (personal  communication  Dr.  Percy  Donaghay).  Despite  these  differences, 
spatial  patterns  of  ocean  color  above  the  sea  surface  always  reflected  similar  modifications 
when  optical  submarine  layering  became  deeper.  In  general,  we  found  in  every  experiment  that 
skewness  of  R 1  switched  from  positive  to  negative,  as  the  optical  submarine  layer  was  placed 
farther  away  from  the  sea  surface  (Fig.  2).  Based  on  two  standard  errors  [i.e.,  ^95  probability  of 
rejecting  null  hypothesis  fi\  =  (i2,  where  /i  is  the  arithmetic  average,  and  subscript  1  and  2  refer 
to  V^l(NS)  and  i/r/?l(FS),  respectively],  we  estimated  the  uncertainty  of  x//R  1  for  each  case 
study  and  found  that  a  minimum  of  15,  5,  and  103  profiles  are  needed  to  differentiate  shallow 
from  “deeper”  layers  in  Black  Sea,  Monterrey  Bay,  and  East  Sound,  respectively.  The  size  of 
these  ideal  datasets  varied  between  1 .3  km  (e.g.,  Monterrey  Bay)  and  9  km  (e.g.,  East  Sound). 
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Fig.  2  Vertical  position  of  subsurface  optical  layers  and  spatial  patterns  of  Rrs.  Depth  variations 
of  a(675)  for  shallow  [(a)— (c)]  and  deep  [(d)— (f)]  case  studies  are  indicated  for  Black  Sea  [(a) 
and  (d)]t  Monterrey  Bay  [(b)  and  (e)],  and  East  Sound  [(c)  and  (f)].  V'Fh  is  a  dimensionless 
variable  defining  the  spatial  skewness  of  flrs(443)/ftrs(490)  (lower-right  corner  in  each  panel), 
two  standard  errors  of  \j/R\  are  indicated  between  parentheses,  FOD  is  the  first  optical  depth 
in  m  (black  dashed  line),  PD  is  the  pycnocline  depth  in  m  (black  solid  line).  Missing  data  are 
depicted  in  white.  Plots  are  based  on  1-m  vertically  averaged  a(675)  data  interpolated  between 
Scanfish  dives. 


Given  the  relatively  shallow  detection  depth  [i.e.,  smaller  than  depth  of  subsurface  a(675) 
layer]  of  passive  sensors  in  these  waters,  the  aforementioned  spatial  changes  in  ocean  color  could 
be  mostly  explained  by  optical  changes  within  the  layer  above  the  main  “opticline.”  Indeed,  we 
found  that  skewness  of  inherent  optical  properties  measured  in  this  layer  and  /?1  may  covary  at 
specific  wavelengths  (Table  1).  For  each  experiment,  deepening  of  the  subsurface  optical  layer 
and  associated  decrease  of  x)/R\  was  always  accompanied  by  a  decrease  on  tf(490)//>(490)  above 
the  main  opticline.  This  pattern  was  not  observed  at  shorter  wavelengths  or  using  single  optical 
variables  instead  of  ratios.  Thus,  we  suggest  that  dissolved  colored  components  (i.e.,  constituents 
strongly  absorbing  in  the  UV-blue  spectral  range)  had  a  smaller  influence  on  x//R\  variability 
with  respect  to  that  associated  with  particulates  (e.g.,  phytoplankton  cells,  detritus,  bacteria). 
Also,  the  above  results  imply  that  \f/R\  variability  is  largely  modulated  by  light  absorption  and 
scattering. 

The  aforementioned  covariability  between  \//R\  and  depth  of  the  optical  layer  was  also 
attributed  to  light  spectral  changes  with  depth.  When  layers  are  relatively  shallow,  the  variance 
in  the  reflectance  at  the  wavelength  that  is  attenuated  faster  dominates  the  variance  in  the 
reflectance  ratio.  When  the  layer  is  relatively  deeper,  this  wavelength  stops  varying  since  light 
is  mostly  attenuated  in  the  layer  nearest  the  sea  surface.  Thus,  the  variance  in  the  reflectance  in 
the  other  wavelength  dominates  the  variance  of  the  ratio.  This  change  exhibits  itself  as  a  change 
of  sign  in  the  skewness  of  R\ .  The  variance  in  reflectance  at  a  given  wavelength  itself  is  driven 
by  changes  in  the  depth  of  the  layer  with  the  added  issue  of  exponential  weighting  of  the  optical 
properties  contributing  to  the  reflectance.  Hence,  a  meter  of  deepening  contributes  differently 
to  variance  than  a  meter  of  shallowing. 

Thermohaline  structure  was  usually  coupled  to  depth  changes  of  pigmented  particulates 
as  inferred  from  a(615)  (Fig.  3).  These  results  emphasize  the  importance  of  stratification  on 
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Table  1  Spatial  statistics  of  opticals  properties  above  the  main  opticline.  For  each  case  study, 
spatial  skewness  is  computed  for  different  lOPs  and  wavelengths  influencing  skewness  of 
F?rs  (440)/ Rrs  (488) . 


Optical  variable 

Wavelength  (nm) 

Experiment 

NS 

FS 

a 

440 

BS 

0.29 

0.11 

MB 

-0.28 

0.69 

ES 

16.8 

1.55 

488 

BS 

1.00 

-0.02 

MB 

-0.91 

0.69 

ES 

1.61 

1.85 

b 

440 

BS 

0.42 

-0.27 

MB 

-0.14 

0.74 

ES 

0.02 

0.11 

488 

BS 

0.29 

-0.20 

MB 

-0.11 

0.79 

ES 

0.19 

2.15 

bla 

440 

BS 

0.64 

-0.32 

MB 

-0.32 

0.26 

ES 

-1.47 

-0.52 

alb 

490 

BS 

2.50 

0.28 

MB 

0.69 

-8  10"5 

ES 

1.98 

0.75 

vertical  distribution  of  phytoplankton,  and  are  consistent  with  other  studies  highlighting  the 
major  association  (i.e.,  >70%  of  cases)  between  the  depth  distribution  of  optical  layers  and  water 
column  stability."  Despite  this  overall  correspondence,  decoupling  between  optical,  biological, 
and  physical  variables  may  occur  (e.g.,  ES  in  Fig.  2)  when  optical  constituents  do  not  behave  as 
passive  tracers  (e.g.,  migrant  phytoplankton),  or  they  form  aggregates  that  escape  the  pycnocline 
barrier.  These  ecological  scenarios  may  introduce  large  uncertainties  when  vertical  localization 
of  submarine  optical  layers  is  based  on  hydrographic  profiles.12  Conversely,  \f/R\  was  very 
sensitive  to  vertical  changes  of  subsurface  optical  layers  for  a  broad  range  of  vertical  mixing 
and  ecological  conditions.  However,  we  recommend  more  comparisons  in  stratified  coastal 
waters  with  different  phytoplankton  communities  before  the  use  of  xf/R\  can  be  generalized  as 
a  noninvasive  optical  proxy  for  screening  depth  changes  on  subsurface  optical  layers. 
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Fig.  3  Biological  versus  physical  mechanisms  modulating  vertical  distribution  of  optical  proper¬ 
ties.  Vertical  variation  (A)  of  a( 675)  and  seawater  density  (p)  per  meter  is  plotted  in  loglO  scale. 
Each  data  point  corresponds  to  the  arithmetic  average  within  1  m  bin.  Ascending  and  descending 
Scanfish  measurements  are  included. 
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